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We discuss combined constraints, coming from the cosmic microwave background shift parameter 
1Z, baryon acoustic oscillations (BAO) distance parameter A, and from the latest type la supernovae 
data, imposed on cosmological models which allow sudden future singularities of pressure. We show 
that due to their weakness such sudden singularities may happen in the very near future and that 
at present they can mimic standard dark energy models. 



PACS numbers: 98.80.Es; 98.80.Cq; 04.20.Dw 

I. INTRODUCTION 

Sudden future singularities of pressure |l| (SFS or type 
II singularities [U, yj) are some exotic-type singularities 
which are allowed to occur in the universe. The in- 
spiration for investigating these new types of singular- 
ity was the observation of high-redshift Type la super- 
novae (SNIa) which over the past decade has provided 
strong evidence that the expansion of the universe is ac- 
celerating [f| and, in particular, the further extension of 
these datasets @, 0] which remain consistent with the 
dark energy existing in the form of phantom Q energy. 
Phantom-driven dark energy leads to a big-rip singu- 
larity (BR or type I according to @) in which all the 
matter is dissociated by the phantom-driven dark energy 
which acts as antigravity in a large and a dense universe 
Q. This behavior is of course different from the stan- 
dard picture of cosmic evolution which allows big-bang 
(BB) or big-crunch (BC) standard types of singularities 
only. Besides, phantom energy violates all the energy 
conditions: the null (gc 2 + p > 0), weak (gc 2 > and 
gc 2 + p> 0), strong (gc 2 + p > and gc 2 + 3p > 0), and 
dominant energy (gc 2 > 0, — gc 2 < p < gc 2 ) (here c is 
the speed of light, g - the mass density in kg m~ 3 , and 
p - the pressure). Apart from big- rip and sudden future 
singularities there are numerous other types such as: gen- 
eralized sudden future singularities (GSFS), finite scale 
factor singularities (FSF or type III) which were tested 
against observations in [Io| . big-separation singularities 
(BS or type IV) and ^-singularities (type V) They 
are characterized by violation of all or some of the energy 
conditions which results in a blow-up of all or some of the 
appropriate physical quantities such as: the scale factor, 
the energy density, the pressure, and the barotropic in- 
dex (for a review see Ref. [3|). In fact, all of these sin- 
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gularities may be investigated by using the higher-order 
characteristics of the expansion of the universe known as 
statefinders [19J. 

Sudden future singularity models result in leaving no 
restrictions on an equation of state p = p(g) of the cosmo- 
logical matter, which allows the unconstrained evolution 
of the energy density and pressure. The nature of a sud- 
den future singularity is different from that of a standard 
big-bang singularity and also from an exotic big-rip sin- 
gularity in that it does not exhibit geodesic incompletness 
and the cosmic evolution may eventually be extended be- 
yond it [HI, EH • The only physical characteristic of these 
singularities is a momentarily infinite peak of the tidal 
forces in the universe. In generalized sudden future sin- 
gularity models this peak may also appear in the deriva- 
tives of the tidal forces It is interesting to note that 
these singularities are in a way similar to yet another 
type, which were termed finite density singularities [lij ]. 
The crucial difference is that finite density singularities 
occur as singularities in space rather than in time, which 
means that even at the present moment of cosmic evolu- 
tion they could exist somewhere in the Universe [15j . We 
will not discuss in detail these finite density singularities 
in this paper since they basically appear in cosmologi- 
cal models without homogeneity. On the other hand, it 
is worth mentioning that the sudden future singularities 
are quite generic, since they may arise in both homoge- 
neous [H and inhomogeneous [3] models of the universe. 

In Ref. [17] we discussed the constraints imposed on 
SFS models which came from SNIa data. We showed 
that those data were consistent with the occurrence of 
a sudden future singularity in the very near future - 
as little as 8.7 million years from now. In Ref. [l8[ 
we extended our investigation to confront SFS models 
with other cosmological data, from the cosmic microwave 
background (CMB) and baryon acoustic oscillations. We 
demonstrated that the class of SFS models introduced in 
[l| was not compatible with current observations, in the 
particular case where the asymptotic behaviour of the 
scale factor close to the BB singularity mimics a dust- 
filled Einstein de Sitter universe. In this paper we extend 
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our investigation further, confronting current cosmolog- 
ical data with SFS models that simulate the behaviour 
of more general flat, barotropic fluid models. The paper 
is organized as follows. In Section [TT| we present sudden 
future singularities in a framework which is appropriate 
to further discussion. In Sections IIIHTVl we discuss cur- 
rent observational data on supernovae, shift parameter 
and baryon acoustic oscillations in the universe. In Sec- 
tion I VII we compare these observational data with our 
SFS models and use them to determine constraints on 
the model parameters and present our conclusions. 

II. THE MODELS 

In order to obtain a sudden future singularity consider 
the simple framework of an Einstein-Friedmann cosmol- 
ogy governed by the standard field equations 

3 /a 2 kc 2 \ 

c 2 / a a 2 kc 2 \ 
»7rG \ a a z a z J 

where the energy-momentum conservation law 

e{t) = -^( e {t) + ^§) , (ii.3) 

is trivially fulfilled due to the Bianchi identity. Here a = 
a(t) is the scale factor, the dot means the derivative with 
respect to time t, G is the gravitational constant, and 
the curvature index k = 0, ±1. What is crucial to obtain 
a sudden future singularity is that in the general case no 
link between the energy density and pressure (i.e. the 
equation of state) is specified (but see below for the case 
of the flat models). 

From equations (|II.1[1 - (|II.2|) one can easily see that a 
pressure singularity p — > oo occurs when the acceleration 
a —> — oo, no matter that the value of the energy density 
g and the scale factor a(t) are regular. Since in that case 
| p |> g, it is clear that the dominant energy condition 
is violated. This condition can be achieved if the scale 
factor takes the form [J 0, 53] 

a(y) = a s [5 + (1 - S) y m -6(1- y) n ] , y = t/t s 

(II.4) 

with the appropriate choice of the constants S, t s ,a s , to, n. 
First publication of the eq. (|II.4[) was in where it was 
used to show that closed universes obeying energy condi- 
tions did not necessarily recollapse. Notice that in equa- 
tion (1II.4[) a s has the unit of length and all the terms 
(including y) in the bracket are dimensionless. before 
the expansion max ie a sudden singularity as I termed 
it later. It is interesting to note that for a flat (k = 0) 
Friedmann model we have an explicit relation between 
the pressure and the energy density, though with a time- 
dependent barotropic index, in the form 

p s (t)=w s (t)g s {t) , (II.5) 



where 

w.(t) = j[2q(t)-l] , (II.6) 

and q(t) = —da/ a 2 is the deceleration parameter. The 
index 's' has been attached to mark the fact that we have 
a type of matter which is of sudden future singularity 
origin and this may play the role of the dark energy. In 
the limit 8—^0 we have q(t) — (I — to) / m which for dust 
(i.e. for to = 2/3) gives q(t) — 1/2, as required. In fact, 
the relations (|II.5[) and (|II.6[) remain of the same form 
if the time t is replaced by the scaled time to a sudden 
singularity y. 

The standard Friedmann limit (i.e. models without an 
SFS) of (|II.4[) is achieved when 5—^0; hence (5 becomes 
the "non-standardicity" parameter of SFS models. Addi- 
tionally, notwithstanding Ref. [l[ and in agreement with 
the field equations 1 [) - 2[) . we assume that 5 can be 
both positive and negative leading to a deceleration or 
an acceleration (cf. pi-4|l ) of the universe, respectively. 

It is important to our discussion that the asymp- 
totic behaviour of the scale factor (111.41) close to the 
BB singularity at t — is given by a simple power- 
law <zbb = y m , simulating the behaviour of flat k = 
barotropic fluid models with m = 2/[3(w + 1)] . This 
allows us to preserve all the standard observed charac- 
teristics of early universe cosmology - such as the cos- 
mic microwave background, density perturbations, nu- 
cleosynthesis etc. - provided we choose an appropri- 
ate value of m. On the other hand, close to an SFS 
the asymptotic behaviour of the scale factor is non- 
standard, asFs(y) = &s [1 — m (1 — <5) (1 — 2/) ], showing 
that asFS = Os for t — t s (i.e. y = 1) at the SFS. No- 
tice that one does not violate the energy conditions if the 
parameter m lies in the range 

0<m<l (w>-l/3), (II.7) 

This range of values is, in fact, equivalent to a standard 
(neither quintessence-like nor phantom-like) evolution of 
the universe. However, with no adverse impact on the 
field equations (|II.1|) - (III.2[) . one could also extend the val- 
ues of to to lie in the complementary ranges [l2| to > 1 
(i.e. —I < w < —1/3) for quintessence, and m < (i.e. 
w < — 1) for phantom models, although these ranges may 
lead to violation of the strong and weak energy condi- 
tions, respectively. 

A special case in which sudden future singularities al- 
low for an explicit equation of state is the anti-Chaplygin 
gas model [2l| 

K*) = "4 G4>0) , (H8) 
g(t) 

and the pressure singularity, known in this case as the so- 
called 'Big-Brake', is achieved for g —> and so p —> oo. 
Here we can easily check that one is able to obtain an 
anti-Chaplygin gas behavior for the scale factor (III. 41) . 
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In order to prove that we need simply to consider the 
first time derivative of (|II.4I) 



a(t) = a s 



m{l-5) m _ x 
1 V 



(II.9) 



We require d — > 0, which corresponds to g — > at y = 1 
so that we have a condition that either m — > or 6 — > 1. 
In fact, these conditions are almost equivalent since 



lim a(y)=a s [l-S(l-y) n ] , 

m— >0 

lima(y) = a s [l-(l-y) n ] . 

o — >•! 



(11.10) 

(11.11) 



It follows that (|II.10j) is more general since additionally 
it does not restrict 51, and moreover it has a standard 
Friedmann limit 5 — > 0, though a static one. It is useful 
to calculate from (|IlT0|) that (cf. pIL4)l ) 



H(y) = 



nS(l -y)™- 1 
[1 - <J(1 - ' 
(n-1) 



(l-2/) r 



(11.12) 
(11.13) 



Using the conservation law 3[) together with the 
barotropic relation 5[) and (|II.6[) (which holds for flat 
fc = models only) we may write 



dg s 
Qs 

which integrates to give 
Qs(t) = ^06-exp 

= gas exp 



-2 [<z(i) 



1] -dt 

a 



2 / [q{t') + l]H(t')dt' 

to 

[w(t') + 1] H{t')dt' 



(11.14) 

(11.15) 
.(11.16) 



The following regimes of the relations (|II.15|) and 
(III.16[) are in order. Firstly, the regime of the early uni- 
verse can be recovered by taking f<f s (which is equiv- 
alent to (5 — > 0) one has 



Qs/gos = (yo/y) 2 oc r 



(11.17) 



Secondly, the near-to-SFS limit can be obtained by tak- 
ing t oc t s which gives 



g s (t) = go s 



(1 - a)t s + at 



(1 - a)t s + at a 



cx g 0s t 2 , (11.18) 



with a = m(l — 5). In fact, one can see from pi.!8[) 
that the contribution from the 'SFS-driven' matter is 
more and more important and starts dominating from 
the current moment of the evolution, i.e., t < t Q implies 
g s < g sQ , t = t implies g s = g Qs , t > t a implies g s > g s0 , 
and t — > t s implies g s = go s t 2 (t s — at s + ato) 



Defining the density parameter 

_8ttG 
6H 



(11.19) 



one can write down the flat Friedmann equation (|II.1[) in 
the form 



a 



a- 3 
H^fl s0 exp 



Q.(t) 



(11.20) 



-3 / [w{t') + l}H{t')dt' 

'to 



H(y) = H Jn s0 exp -3/ [w(y>) + 1] H{y>)dy> . 

V L Jyo 

(11.21) 

In this paper we consider three observational con- 
straints on SFS models: luminosity distance moduli to 
SNIa, the CMB shift parameter, which is a scaled dis- 
tance to the last scattering surface of the cosmic mi- 
crowave background, and the distance parameter con- 
strained by baryon acoustic oscillations. 



III. SUPERNOVAE 

We proceed within the framework of Friedmann cos- 
mology, and consider an observer located at r = at 
coordinate time t = to- The observer receives a light 
ray emitted at r = r\ at coordinate time t — ti and, 
according to (|II.4|) . with redshift given by 

i , afro) _ 5+(l-6)yF-6(l-y Q ) n 
a(h) 5+(l-6)y?-6(l- yi ) n ' 

where yo = y(to) and y\ = y(ti). We then have a stan- 
dard null geodesic equation 



dr 



o VI - fcr 2 J tl o(t) ' 



to 



cdt 



(III.2) 



with the scale factor a(t) given by (|II.4j) . For a flat Fried- 
mann model we can write down the radial coordinate to 
an observer in any of the forms 



ri 



*° cdt 
cdz 



Vo (ct s )dy 

vi a (y) 



cda 



cda 



'"' J ai Ha 2 



H(z)a H a J E(z) 



dz 



(III.3) 



where we have used that y — t/t s and so dy = dt/t s 
which implies a'(y) = da(y)/dy — t s da{t)/dt = t s a(t), 
"■'(yo) = t s d(t ) and 

H = H(t ) = = —H(yo) ; H(y) = — — . 

t s ao ts a[y) 

(HI.4) 

In (|III.3[) the transition from the integral of da to the 
integral of dz was given by the application of the defini- 
tion of redshift 1|) . Besides, due a lack of an analytic 
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form for the equation of state for SFS models the func- 
tion E(z) can only be given by a formula which involves 
an integral over z, as follows 



E\z) 



H 2 (z) 
Hi 



Sl<5o(l + z) 3 exp 



dz 



,2«(*') 



1 



1 



(III.5) 

in (|III.3|) with Q$o being the density parameter of sudden- 
future- singularity- driven dark energy [3]. It is easy to 
notice that in the limit S — > and m — 2/3 one has 
that E(z) = r2ao(l + z ) 3 i as f° r the standard matter 
dominated case. Let us recall that the standard formula 
for the models which include the dark energy component 
£l w o reads as H^, [2(| 

e 2 (z) = n m0 (i + zf + n w0 (i + z) 3 (- +1 ) + n M (i + zf , 

(III.6) 

where Qko is the curvature component and p — wg. How- 
ever, in our further calculations we will not be expressing 
E(z) in either forms ()III.5|) or (|III.6|) . using an explicit 
form of a(y) as in (|III.3|) instead. 

The luminosity distance to e.g. a supernova observed 
at redshift, z, is given by 



d L (z) = (1 + z)a(t )n 



(III.7) 



where r\ is given in one of the forms (|III.3[) . For our 
calculations we have applied the following expression 



d L (z) = (1 + z)a(y )ct s 



yo _dy_ 

>vi a (y) 

The distance modulus is 

fi(z) = 51og 10 d L (z) + 25. 



(III.8) 



(III.9) 



We compared the luminosity distance predicted for SFS 
models with SNIa data using the SCP Union2 
dataset, consisting of 557 supernovae, which is the largest 
compilation published to date. 



IV. SHIFT PARAMETER 

The standard formula for the CMB shift parameter is 
given by [HHi]: 



n 



ITT 



jTT 

n 



(IV.l) 



where lf T is the temperature perturbation CMB spec- 
trum multipole of the first acoustic peak in the model 
under consideration and l'^ T corresponds to a reference 
flat standard Cold Dark Matter (CDM) model. The mul- 
tipole number is related to an angular scale of the sound 
horizon r s at decoupling by [2(| 127) 



^- cx - 

dA h 



(IV.2) 



where 



r s = adccS(r s ) = adccS yj c. 



a{t) J ' 



(IV.3) 



with c s being the sound velocity and the angular diameter 
distance reads as 



f to dt 
/tdcc a (0 



(IV.4) 



where c is the velocity of light and S[r) = r for k = 0. 

Following 25] we can write, using IIV.1I and IIV.2I as 
follows 



K = 



r s d' A (z' dcc ) 
r' s d A {z dcc ) 



(IV.5) 



which, by assuming that at decoupling the amount of 
radiation was the same in both the flat reference standard 
CDM model and in our SFS model (which we assume to 
be just the same as a standard matter-radiation model of 
the early universe, since SFS models do not change the 
evolution there) we have that 



1 



(IV.6) 



r 's V"mO 

On the other hand, for a reference standard CDM model 



d , = 2ca' dc 



a H a 



«o 



2ca' 



dec 



a Ho 



f(O. r0 , a' dec ) 



(IV.7) 



while for our SFS model the angular diameter distance 
is given by 



d-A = Odec^de 



(IV.8) 



with rd c given by (|IIL3|) taken at decoupling. Using the 
above, we may write that for our SFS models the shift 
parameter is 



n 



2c 



2cf. 



Hoao^/Qmordec a'(yo)v / ^mO?'do 



(IV.9) 



where we have assumed that the function f(£l' r0 , a' dec ) is 
approximately unity pEl ]. 

Finally, the rescaled shift parameter is 



n 



V"roO r de 



\/n m0 a'(y ) 



dy 



dz 

W) 



a (y) 



(IV.10) 



where in order to obtain the last expression we have used 
dnL4| . The WMAP data gives K = 1.70±0.03 Note 
that, as we discussed in [l^l, this value is not directly 
"observed" but is derived from the CMB data assuming 
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a specific class of cosmological model; hence one must 
be careful in employing it to test our SFS model. How- 
ever, as we showed in [18|], the effective equation of state 
for our SFS model displays general similarity to that of 
the concordance model and on this basis we consider the 
use of the "observed" value of the shift parameter to be 
appropriate for our SFS model too. 

Moreover in [l8[ we also considered the acoustic scale, 
l a , and how its incorporation into the analysis could re- 
sult in parameter constraints which gave a closer approxi- 
mation to those obtained from a fit to the full CMB power 
spectrum. However, upon varying the parameter m from 
the Einstein-de Sitter value of 2/3 considered in [181 ] . it is 
not clear that the use of the acoustic scale remains valid. 
To see this we recall that the acoustic scale is defined in 
terms of two distinct distance scales: it is proportional 
to the ratio of the angular diameter distance, cLa, and 
the sound horizon, r s , both evaluated at the redshift of 
recombination, zqmb- 



_<Ia{zcmb) 



r s {zcMB) 

where the comoving sound horizon is defined as: 
cAz' 



E(z') 



with 



1 



3(1 +R b a) 



(IV.ll) 



gV.12) 



where E(z) = H(z)/Hq and c s is the speed of sound 
with R b = 31500n b h 2 (^f-)- 4 ,T CM B = 2.7K. Now the 
potential difficulty in carrying out the l a calculation for 
the m 7^ 2/3 case arises from the sound speed calcula- 
tion. Previously in the to = 2/3 case we could interpret 
the sound speed as being set at the time of recombina- 
tion in the "standard early universe" . The value of the 
acoustic scale would still in general be different for our 
SFS universe due to the SFS model's different expansion 
history and hence different angular diameter distance- 
redshift relation. In the case of m / 2/3, however, the 
above simplification is no longer valid and both the sound 
horizon and angular diameter distance will in general be 
different for our SFS universe. In particular computing 
the sound speed in the m ^ 2/3 case would require more 
detailed modelling of the SFS fluid in the early universe 
which lies beyond the scope of this paper. Hence we do 
not consider the acoustic scale further in this paper. 



V. BARYON ACOUSTIC OSCILLATIONS 

The Alcock-Paczyhski effect [29| states that one is able 
to calculate the distortion of a spherical object in the sky 
without knowing its true size. This can be done by mea- 
suring its transverse extent using the angular diameter 
distance, r 



I 

A9 ' 



(V.l) 



where I and AO are the linear and angular size of an 
object, and its line-of-sight extent, Ar, using the redshift 
distance 



Ar = 



cAt ct s Ay 
a{t) a(y) 



(V.2) 



(see e.g. Ref. [25j). As a result one can define the volume 
distance, Dy, as 



Dy = r 2 Ar 



so that using (|III.3j) . one has 



Vo ct s dy 

vi a (v) 

c f z dz 



ct s Ay 

a(y) 



a H Jo E(z)J \a H E(z 



Az 



(V.3) 



(V.4) 



(V.5) 



Ref. [30jj uses the acoustic peak scale for 46748 luminous 
red galaxies selected from the Sloan Digital Sky Survey 
to measure Dy(Az = zbao = 0.35) = 1370 ± 64 Mpc. 

Usually, it is more convenient to work with a dimen- 
sionless quantity A which for our SFS model (|II.4[) is 
obtained multiplying (IV.5I) by \Z£l m o(aoHo) / (czbao) to 
obtain 



A = n 1 i 2 l E( ZBAO ) 



-1/3 



1 


f Zl dz ' 


.£bao J 


o E(z)\ 



2/3 



(V.6) 



or, using the definition of E(z) given by the first equality 
in (|III.5j) . together with pil.4j> . 



A = 



\/fL m oa'(y ) 

[o'(?/BAo)]3 



a(yo) 



1 



dy 

ZBAO Jj/bao 

a(y) 



(V.7) 



The same result can be obtained from (IV.4I) if we com- 
pute the difference Ay = yo — Vbao and a(y) = a(yo) in 
(ED, i.e., 



Dv = ct s 



yo - Vbao 



a(yo) 



Vo dy 
y B Ao <V) 



(V.8) 



and use the fact that ct s = c/H a'(yo)/a(yo)- Following 
Ref. [3(| the parameter A should have the value 

(ji \ —0.35 
09s) ±0 -° 17 ' (V - 9) 

where n is the spectral index (now taken to be ~ 0.96). 

VI. RESULTS AND CONCLUSIONS 

We used a Bayesian framework to confront our SFS 
model with the cosmological observations discussed in 



6 



the previous sections. For each cosmological probe we 
took the likelihood function to be Gaussian in form, i.e. 



1 



p(data|0) cx exp(— ~x ) 



(VI.l) 



where denotes the parameters of the SFS model and 
"data" denotes generically the observed data for one of 
the three cosmological probes. For the SNIa data % 2 
takes the form 



2 _ {^obs( z i) ~~ Hprcd( z i)) 2 

xsn — 

i=l 



a 



(VI.2) 



ml 



where o~i is the quoted observational error on the i th 
Union2 SNIa and a ln t is the SNIa intrinsic scatter. Fol- 
lowing [lU we take cr in t = 0.15. For the CMB shift pa- 
rameter \ 2 takes the form 



2 _ (ft -1.70)' 

Xn ~ 



0.03 2 



(VI.3) 



while for the BAO distance parameter it takes the form 



x\ 



(A- 0.469) 2 
0.017 2 



(VI.4) 



Since our three cosmological probes are mutually inde- 
pendent, their joint likelihood function is given by the 
product of their individual likelihoods, and thus takes 
the form 

p(all data|8) oc exp(— ^Xtot): (VI-5) 

where Xtot = Xsn + Xn + Xa- 

We used Bayes' theorem and a Markov Chain Monte 
Carlo (MCMC) approach [20| to obtain posterior prob- 
ability distributions for the SFS model parameters: n, 
to, S and ya = to/t s , where to is the present age of the 
universe. We employed the Metropolis-Hastings algo- 
rithm, adopting uniform priors for each of the param- 
eters: S € (-30, 0), n e (1, 2), to G (0, 3), y € (0, 1). 

In Fig. Q] we present contour plots showing the joint 
marginal posterior distribution for each pair of SFS 
model parameters. Each sub-panel shows three contours, 
denoting 68.3%, 95.4% and 99.73% (from light gray to 
dark grey respectively) credible regions. 

In Fig. [5]we also show marginal posterior distributions 
for each of the model parameters individually, again de- 
termined via Markov Chain Monte Carlo. It is immedi- 
ately clear from Figs. [T] and [5] that the current cosmologi- 
cal data strongly exclude our SFS model for to = 2/3 (in- 
dicated by the dashed horizontal and vertical lines in the 
relevant sub-panels of each figure) . In the leftmost panel 
of Fig. [2J for example, we see that the marginal pos- 
terior distribution for to is strongly peaked between 0.7 
and 1.0, but is negligible outside the range (0.7, 1.0). In 
other words our SFS model is incompatible with current 
observations in the particular case where the asymptotic 
behaviour of the scale factor close to the BB singularity 



mimics a dust-filled Einstein de Sitter universe. This of 
course is not surprising since as noted in Section U it was 
the conclusion we reached in [l8j]. On the other hand 
it is also clear from Figs. [T] and [5] that there is a sig- 
nificant region of parameter space where the marginal 
posterior distributions of the SFS model parameters ap- 
pear to overlap with each other - i.e. one can identify 
sets of model parameters which occupy the 68.3% credi- 
ble region for all of the marginal posterior distributions 
shown. This does not in itself guarantee that these sets of 
model parameters give an acceptable fit to the SFS model 
since the marginal posterior distributions are projections 
of the fully, four-dimensional joint posterior: the appar- 
ent overlap of the two-dimensional marginal distributions 
for each pair of parameters, for example, could in princi- 
ple be merely a projection effect. We have checked this 
rigorously, however, via the following approach. We con- 
sidered each of the three cosmological probes separately 
and produced MCMC chains that explored the full four- 
dimensional SFS model parameter space for each probe. 
We then identified the MCMC chain points that sampled 
e.g. the 68.3% credible region for the SNIa chain and 
verified (via direct computation of the likelihood values) 
that there was significant overlap in four dimensions be- 
tween these points and the corresponding credible regions 
of both other probes. Thus we established that the ap- 
parent overlap of the marginal posterior distributions for 
the joint likelihood of all three cosmological probes was 
not merely a projection effect. 

It is clear from this analysis, therefore, that the 
SFS model under consideration is compatible with cur- 
rent cosmological data, based on the three cosmological 
probes that we considered. In particular, provided that 
to > 0.72 we see that we can identify ranges for the 
other parameters for which the predicted values of the 
cosmological probes are in good agreement with current 
observations. As the value of m increases towards unity 
the credible regions for the parameters 6 and n extend 
to fill their allowed prior ranges, while the credible re- 
gion for j/o is pushed towards lower values - i.e. the SFS 
model is still compatible with current data but the epoch 
of the sudden future singularity is pushed further into the 
future. Interestingly, however, we see that as m — ¥ 0.72 
the credible regions for n and 5 become much narrower, 
constrained to lie ever closer to n = 2 and 6 = 0, while 
the credible region for yo tends towards unity, so that the 
SFS may happen in the very near future. 

Thus we conclude that our SFS models are compat- 
ible with current observations provided the parameter 
to, which characterizes the near-to-big-bang evolution 
of the scale factor, is at least 0.72. This excludes the 
Einstein de Sitter dust solution, and requires a form 
of matter which has slightly negative pressure (since 
to = 2/3(w + 1), so that w < -0.083). Note also that, 
in order to match the current observations, the value of 
the "non-standardicity" parameter S is seen necessarily 
to be negative. This naturally eliminates the validity of 
anti-Chaplygin gas models [21| given in a special form 
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FIG. 1. Marginalized contours for pairs of parameters are plotted. There are three levels of credible region shown: 68.3%, 
95.4%, 99.73% (from light grey to dark grey respectively) calculated from A, 1Z, and SN la jointly 
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FIG. 2. Unnormalized probability densities for all parameters, calculated for A, 1Z, and SN la jointly. 



by ([II. Ill) which require 8 — > 1 (but it does not elimi- 
nate the models in the form of 10[) despite their static 
Friedmann limit and also other types of anti-Chaplygin 
gas models studied in Ref. |22j). 



In conclusion, similarly as in our previous paper (l7| 
where we applied supernovae data only, here we have 
shown that a sudden future singularity may happen in 
the near future of the universe. 
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